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ABSTRACT 

We present high-resolution X-ray spectroscopy of GRB 020405 obtained with the 
Low Energy Transmission Grating Spectrometer (LETGS) on board the Chan- 
dra X-ray Observatory starting 1.68 days after the burst. The spectrum appears 
featureless, with no evidence for emission lines, absorption edges, or narrow ra- 
diative recombination continua. The continuum can be fitted by a power law of 
photon index F = 1.72 ± 0.21 and temporal decay index a = 1.87 ± 0.1, with a 
marginally significant excess column density of cold gas Nh={4:.7 ± 3.7) x 10^^ 
cm~^ at the redshift of the host galaxy. The absence of iron lines indicates 
that the density of nearby surrounding material was unlikely to be very dense 
(n ^ 5 X 10^^ cm~^) at the time of the Chandra observation. In the case of 
recombination following photoionization in an optically thin medium, most ionic 
species would be completely stripped at lower gas densities than this. In the case 
of a power-law spectrum reflecting off a "cold" , opaque medium of low density, 
negligible emission features would be produced. Alternative to these possible ex- 
planations for the lack of emission features, any X-ray line emission taking place 
in a dense medium in a "nearby reprocessor" scenario might have been over- 
whelmed by the bright afterglow continuum. Although the absence of discrete 
features does not unambiguously test for a connection between GRB 020405 and 
nucleosynthesis, it emphasizes the need for high-resolution X-ray spectroscopy 
to determine the exact emission mechanism responsible for the reported discrete 
lines in other GRB afterglows. 

Subject headings: gamma rays: bursts — X-Rays: general 



1. Introduction 

The nature of gamma-ray burst (GRB) progenitors remains a mystery despite the lo- 
calization of ~ 40 X-ray afterglows^. An increasing number of X-ray observations reporting 
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discrete spectral features in GRB afterglows (Piro et al. 1999; Antonelli et al. 2000; Piro 
et al. 2000) has motivated attempts to find an explanation in terms of GRB progenitors. 
Thus far the claimed features are consistent with n = 2 ^ 1 transitions in all possible charge 
states of Fe, i.e., from Fe K-shell fluorescence at 6.4 keV to H-like Fe Ly a at 6.95 keV. 
More recently Reeves ct al. (2002) have reported the detection of H-like emission from mul- 
tiple a-elements including Si and S in the spectrum of GRB 011211; however, the statistical 
significance of those features is debatable (Rutledge & Sako 2002). The possible presence 
of Fe and/or a-elements in the X-ray spectra of GRB afterglows could provide valuable 
information about the physical conditions in the vicinity of the progenitors. 

Determining the relevant X-ray reprocessing mechanism in GRBs will involve the physi- 
cal properties and geometry of the surrounding medium (see for instance Lazzati, Campana, 
& GhiseUini 1999; Paerels et al. 2000). Discrete spectral features are expected from hot 
gas when electron collisions cause excitations from the ground state. Emission lines can 
also arise via recombination if continuum photons (either from the burst itself or its after- 
glow) create a highly photoionized plasma. In certain cases fluorescence may be important 
if ions retain sufficient numbers of bound electrons. An important diagnostic that might 
discriminate between the different alternatives is the appearance of radiative recombination 
continua (RRC) produced by the return of unbound electrons to the ground state. RRC 
from a photoionized plasma resemble narrow emission lines, and could be used to verify that 
recombination is occurring. A 10 ks grating observation of GRB 991216 (Piro et al. 2000) 
showed a possible iron RRC that could represent the first evidence of a recombining plasma 
near a GRB. Nevertheless, and despite these efforts, the X-ray mechanism responsible for 
discrete features in GRBs is still uncertain due to limited spectral resolution and sensitiv- 
ity of the existing observations (Paerels et al. 2000). In this letter we present results of 
the LETGS/ACIS-S spectroscopy of GRB 020405, and we show how density and geometry 
regulate the feasibility of detecting discrete lines in X-ray afterglows. 



2. X-ray Observation of GRB 020405 

GRB 020405 triggered the Third Interplanetary Network (IPN) on UT 2002 April 5.0288 
with a 25-100 keV fluence 3 x 10~^ erg cm^^ (Hurley et al. 2002) making it at the time the 
highest fluence above 25 keV for a localized burst since GRB 010629 (Ricker et al. 2001). The 
localization of a bright optical transient (Price, Schmidt, & Axelrod 2002a) and a relatively 
low foreground Galactic column density of Ah = 4.3 x 10^° cm~^ along this line of sight 
generated immediate interest in obtaining high-resolution X-ray grating spectroscopy of its 
afterglow. 
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After a timely response by the Chandra X-ray Observatory, the observation of GRB 
020405 started on UT 2002 April 6.711 using the LETGS in conjunction with the ACIS 
detector for a total usable exposure time of 50.6 ks before the satellite went into Earth 
occultation near perigee. At the time of the final Chandra trigger we still lacked knowledge 
of the redshift, therefore the observation was designed with a displacement along the grating 
dispersion axis, Y-OFFSET=2.'45. The offset placed the Oth-order image on the front- 
illuminated ACIS-S2 chip (Figure 1) while the ACIS-S3 and ACIS-Sl backside-illuminated 
chips, with higher quantum efficiency, covered higher spectral orders in the wavelength range 
1.7 A < A < 28.9 A and 28.5 A < A < 55.8 A respectively. A priori, this configuration includes 
redshifted H-hke hnes of Ni, Co, Fe, Ca, Ar, S, Si, Al, Mg, Ne, and O as well as a few L-edges 
and He-hke lines from a number of metals assuming an average GRB redshift^ z ~ 1.3. 



3. Temporal and Spectral Analysis 

The standard Chandra data products were processed with the destreak tool to correct 
a minor flaw in the readout of the chips^. Subsequently, the instrument and mirror response 
were generated using the "makearf" and "makermf" tools. The final 0th and Ist-order 
spectra were then extracted using the proper grating effective area function (ARF) and 
redistribution matrix file (RMF) for each setting. For the timing analysis we used counts 
extracted from a circular aperture centered on the Oth-order image (we obtain consistent 
results using the dispersed data). A light curve of the LETGS events as a function of 
time shows a decay of the X-ray afterglow of GRB 020405 (Figure 2). The flux can be 
characterized as a power law oc t"'^ (where t stands for time after the initial burst event) 
with a = 1.87 ±0.1. 

The spectral analysis was undertaken using the spectral fitting package XSPEC. We 
analyzed the -|-1 and -1 orders individually as well as coadded. The results obtained in 
each case are statistically consistent. A fit of the coadded Ist-order extracted spectrum is 
well represented by an absorbed power law with photon index F = 1.74 ± 0.22 and column 
density A^h = (1-94 ± 1.14) x 10^^ cm'^ {xl = 0.26). Rather than leaving A^h as a free 
parameter, we also fitted a power-law model with Galactic and redshifted absorption at the 
host redshift, z = 0.690 (Masetti et al. 2002; Price et al. 2002b) by treating the Galactic 
absorption along the line of sight. Ah = 4.3x10^° cm^^, as a fixed parameter. The best 
fit for such a model has F = 1.72 ± 0.21 with an absorbing column density at the host 
A^H(host) = (4.7 ± 3.7) x lO^^cm"^ {xl = 0.26). Figure 3 shows the coadded Ist-order 
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spectrum and power-law fit including residuals. The X-ray flux Fa;(0.2-10 keV) = (1.36 
± 0.25) xlO~^^ ergs cm~^, measured 1.71 days after the burst, with a normalization 
constant A = 2.41 xlO~'^ photons keV^^ cm~^ s~^ at 1 keV is similar to the fiux reported 
for the bright GRB 010222 at the same epoch (in't Zand et al. 2001). Although the value 
of A^H is not highly constraining, it suggests that absorption at the host is non- negligible. 

The X-ray fiux, in comparison with optical photometry obtained around the same epoch, 
1.7 days after the burst (Bersier et al. 2002, Price et al 2002b), can be used to derive a 
broad-band slope f3ox = 0.74. Remarkably, this is identical to the X-ray spectral index itself, 
f3j. = r — 1 = 0.72 ±0.21. However, a smooth extrapolation through the BVRI photometric 
points is not indicated, since fits to the BVRI spectrum yield slope Po = 1-25 (Price et 
al. 2002b) or = 1.43 (Bersier et al. 2002), the later even after correcting for Galactic 
extinction. Thus, the optical slope is steeper than the X-ray, seeming to require somewhere 
a concave upward inflection in the spectrum, which is inconsistent with basic synchrotron 
afterglow theory. This kind of discrepancy is common in afterglow spectra (e.g., Halpern et 
al. 1998), and is most simply understood as requiring additional dereddening of the optical 
spectrum to account for local extinction in the host. In this case, the marginally significant 
A^H(host) derived from the fit to the X-ray spectrum, if accompanied by typical amounts of 
optical extinction, would be enough to restore a concave downward shape to the broad-band 
spectrum. 

Alternatively, this broad-band spectrum can be described as having an X-ray excess, 
which would require a mechanism in addition to synchrotron emission to be operating. For 
example, an X-ray contribution from inverse Compton scattering may be indicated, as was 
deduced by modeUing of the broad-band afterglow of GRB 000926 by Harrison et al. (2001). 

Of the simple synchrotron models involving expansion into a uniform-density medium 
(Sari, Piran, & Halpern 1999), the observed combinations of spectral slope and temporal 
decay in the optical and X-ray are most consistent with the evolution of a jet-like afterglow 
after both the cooling break and the jet break have passed. However, the required electron 
energy power- law slope p would be slightly less than 2, which violates the assumptions of the 
model from which the relations a — p and P — p/2 were derived. A different environment in 
which the medium has a stellar-wind density profile n oc should be considered, because 
its predictions allow a more reasonable value of p. In the adiabatic case where the cooling 
frequency has not yet passed, a = {3p — l)/4 = {3(3 + l)/2 (Chevalier & Li 2000). For 
example, an acceptable fit in the wind scenario has p = 2.7, a = 1.78, and (3 = 0.85. A wind 
model should be of interest in the context of evidence for a massive progenitor star from the 
possible detection of a supernova in the late-time decay curve (Price et al. 2002b). 

The Ist-order spectrum for the whole observation shows no evidence for discrete emission 
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features, absorption edges or narrow RRC signatures across the spectrum. The edge seen 
around 5.62 A corresponds to an instrumental absorption edge in gold material used in the 
LETGS. In addition to searching for features in the summed observation, the data were 
divided into equal time intervals of various durations (AT= 5.1 ks, 10.1 ks, 16.87 ks) to 
test for transient features that might have been present for part of the observation, as was 
reported in the case of GRB 970508 (Piro et al. 1999). The search for transient features did 
not produce any detection at the 3 a level. We also looked for evidence of discrete features 
in the Oth-order spectrum. Although there are no features evident, the detection limits are 
not as sensitive as in the dispersed spectra because of the reduced spectral resolution and 
efficiency of the ACIS-S2 chip. 

In order to quantify the absence of discrete features, we make a simplified assumption 
and adopt a fixed discrete line width (Trest = 0.46 keV in the GRB rest frame (iTobs = 0.27 
keV z — 0.690). This value is borrowed from the Fe hue measurement in the GRB 991216 

spectrum as reported by Piro et al. 2000. We then proceeded to fit Gaussian line profiles 
at their redshifted energies for multiple transitions and used the power-law continuum level 
to derive the line upper limits listed in Table 1. Finally, we determined upper limits for the 
rest-frame equivalent widths at the redshifted location of the transitions listed in Table 1. 

4. X-ray Line Emission Mechanisms 

The main objective of this grating spectroscopy was to acquire a sensitive soft X-ray 
spectrum to derive constraints on material in the proximity of the GRB. This is in fact the 
brightest case where X-ray spectroscopy has been obtained for a GRB in which a SN-like 
signature in the late-time decay curve (Price et al. 2002b) has been reported. A significant 
detection of Fe peak elements and particularly medium-Z a elements (C through Si), where 
the highest spectral sensitivity for this choice of instrument is achieved, could have helped 
estabhsh the nature of GRB progenitors. Unfortunately, the absence of discrete features in 
the spectrum of GRB 020405 does not provide a direct hnk between GRB 020405 and nucle- 
osynthesis. Nevertheless we use the observed X-ray spectrum to place tentative constraints 
on the relevant emission mechanisms. 



4.1. Collisional Excitation 

There has been a resurgence of interest in the possible role of coUisional excitation as a 
line emission mechanism in GRBs. The reason is mainly the MEKAL plasma model fit to 
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the reported spectral emission lines from GRB 011211 (Reeves et al. 2002). A brief review 
of the thermal history of GRB 020405 is needed to evaluate the importance of collisional 
excitation to its X-ray afterglow emission. During the first 60 seconds, the burst bathes 
the progenitor gas with photons of average energy E ^ 0.1 MeV (Price et al. 2002b). For 
the prompt spectrum and peak energy of the burst Ep (Band et al. 1993), the photon 
temperature averages Tprompt ~ 10^ K. However, we note that the Compton time scale 
required to equilibrate is much longer than the duration of the burst unless the electron 
number density n 10^^ cm~^ (Paerels et al. 2000). In other words, unless the medium is 
very dense, the Compton temperature of the prompt emission will not be reached. Current 
numerical models of proposed progenitor collapse have considered central densities of the 
required order (MacFadyen, Woosley, & Heger 2001). Nonetheless, even if occurring at early 
times, such a transient effect will probably be short lived because the plasma cools as the 
afterglow begins. In fact, the plasma temperature will settle to a new value corresponding to 
the balance between Compton heating and cooling and determined by the afterglow spectral 
parameters. Consequently, in the nonrelativistic limit, the temperature of the emitting 
material should not exceed the afterglow Compton temperature Tq given by Ross, Fabian, 
& Young (1999) and expressed as 



where ue oc E^-'^, E^ = 0.1 keV, and E2 = 100.0 keV gives Tc - 4x WK for F = 1.72. 

Additional losses due to bremsstrahlung and atomic transitions continue to cool the 
plasma below Tc as the cooling time scale becomes comparable to the duration of the af- 
terglow. The absence of statistically significant emission in the X-ray spectrum implies that 
collisional excitation does not play a significant role in driving line emission by the time of the 
X-ray observation. In short, collisional excitation by afterglow photons might be important 
only when there is a sufficiently enough dense medium (n ^ 10^° cm^) for Compton inter- 
actions to equilibrate, and if an additional heating mechanism is present that can maintain 
the gas around temperatures T ^ 10^ K, as suggested by Paerels et al. (2000). 



As the ionizing afterglow continuum encounters the external medium, discrete emission 
might also be excited by recombination or fluorescence. Depending on the optical depth 
and geometry of the medium, we can refer to this mechanism as either transmission or 
reflection (Vietri et al. 2001). In transmission, an optically thin medium (r < 1) might 
produce photon cascades following recombination. In reflection, emission can be produced 




(1) 



4.2. 



Recombination &; Fluorescence 
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by the afterglow spectrum reflecting off an optically thick slab (r S> 1). Although the hne 
equivalent widths of order a few keV reported in prior observations favor reflection over 
transmission, we will consider both possibilities here. The structure of recombination in 
either case will be determined by the ionization parameter C, = L^/nR'^ (Tarter, Tucker, 
& Salpeter 1969) where is the ionizing luminosity, n is the number density and R the 
distance from the ionizing source. 

To start, we derive the constraint on density for recombination taking place in an op- 
tically thin, homogeneous medium. Nominally, since the Chandra observation ended 2.27 
days after the GRB any emitting material must have been located within a radius R of the 
GRB where 

^ ctobs 1^ ^ ^ 1^ cm, (2) 

{1 + z) (1 -cos^) (1 -cos^) ' ^ ' 

where 6 is the angle between the observer and the emitting material as seen from the GRB. 
Furthermore the light-distance traveled in the rest frame in the duration of the X-ray obser- 
vation is AR = 8.98 X 10^"^ cm. 

Next we adopt a time-dependent X-ray luminosity 



L^(0.2 - 10 keV) 3.3 x 10^^— ^— ergs s^^ (3) 

^ ' An Vl-71daysy ^ ^ ^ 

with energy index f3x = 0.72 where Q/Att corresponds to the jet solid angle. Price et al. 
(2002b) estimated a jet opening angle of (5.°83 ± 0.°69) n^/^ for GRB 020405 where n is the 
particle density in the nearby medium. However, since the solid angle of the X-ray continuum 
is poorly constrained at the time of our observation, which may overlap the "jet break," we 
leave Q as a free parameter in the following discussion. At z = 0.690, GRB 020405 has a 
luminosity distance — 4505 Mpc, assuming a. Hq — 65 km s~^ Mpc~^, = 0.3, D,\ — 0.7 
cosmology. Hence 

^ 3.3 X 10^ Q / t 

f ~ 7^ 

iT-ufie ^'^ V 1.71 days 

where rig is the distance from the ionizing source in units of 10^^ cm and rin the number 
density in units of 10^^ cm~^. 

As noted above, the state of the gas, including the ionization fraction / for each species, 
is largely determined by the ionization parameter ^. In order to estimate / more accurately 
as a function of ^, we ran the XSTAR code (Kallman & Bautista 2001). XSTAR produces 
ionization fractions for H, He, C, N, O, Ne, Mg, Si, S, Ar, Ca, Fe, and Ni using the appropriate 
photoionization and recombination rates. For our model we adopted an input with power-law 
energy index /3j.=0.72 and a homogeneous medium with solar chemical composition within 
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the radius of emission R. The output of XSTAR shows that the ionization fraction of bound 
Fe"*"^^ starts to be significant near ^ ~ 10^, reaching a maximum at ^ ~ 10^. Other species 
require lower ionization parameters for significant ionization stages with bound electrons to 
be present. 

In order to constrain the density of the medium, we express the recombination luminosity 
of any line L^ne as 

^ ~ innARa{Z, T)fAz^ (5) 



where Az is the element abundance, E'lme is the line energy, / the ionization fraction, a{Z,T) 
the recombination coefficient (Seaton 1959). Assuming ^ < 10^, an ionization fraction of 
order 0.1 (possibly overestimated), and standard solar abundances, the upper limit to the 
Fe Ka line luminosity listed in Table 1 translates into an upper limit in density 

n^5 X lO^^l^ cm-3 (6) 

This expression shows that significant ionization fractions require densities n ^ 5 x 10^^ cm~^ 
for species to be recombining under conditions of ionization equilibrium. For lower densities, 
most species including Fe would be completely stripped, which would account for the absence 
of discrete features in the spectrum of GRB 020405. Moreover, the medium becomes optically 
thick near n ^ 2 x 10^ cm~^. The range of ionization fraction maxima for Fe"^^^ and other 
species is relatively narrow and could serve as an indicator of the recombination mechanism 
in transmission as ^ evolves with time due to the afterglow decay. Another caveat in this 
analysis is our assumed standard solar abundances. Supersolar compositions (10-100 times 
solar) would relax the restrictions on density by the appropriate factor, as the line luminosity 
is linearly dependent on both abundance and density. 

The reflection model for an afterglow spectrum off a "cold", optically thick medium 
has been investigated numerically by a number of groups (Ballantyne & Ramirez- Ruiz 2001; 
Lazzati, Ramirez- Ruiz, & Rees 2002). In particular, the reflection model has been put forth 
to explain the X-ray spectrum of GRB 991216 (see for instance Vietri et al. 2001). Using 
XSTAR once more with the estimated power law spectrum, solar abundance and radius of 
emission R, we find that refiection models appear most efficient for ionization parameters 
^ ^ 10^. This result is consistent with previous work (Ballantyne & Ramirez- Ruiz 2001; Kail- 
man, Mcszaros, & Rees 2001). When fine emission from Fe and other elements is absent, any 
geometry with the optimal ionization parameter (^ ^ 10^) implies a particle density of order 
n ^ 2.7 X 10^° cm~^ at the assumed emitting region R. Lower column densities increase the 
ionization parameter and the quoted models tend to become near-perfect reflectors thus ex- 
hibiting negligible emission features (Ross, Fabian, & Brandt 1996). A near-perfect reflector 
is also consistent with the absence of emission features in the spectrum of GRB 020405. 
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Below the ionization stage Fe+^^, Fe Ka fluorescence emission might follow photoion- 
ization. Although the fluorescent line luminosity is typically independent of the ionization 
parameter ^, fluorescence is usually important when ^ < 10^ (Kallman, Meszaros, & Rees 
2001). This implies that fluorescence might only occur if n ^ 2.7 x 10^^ cm^"^ for emission 
material within R. The non-detection of Fe Ka at 6.4 keV appears to rule out significant 
fluorescence in a solar abundant medium. 



4.3. Afterglow Continuum: The Photon Curtain Effect 

Two general classes of models have been put forward to explain the reported GRB 
X-ray features. In the flrst instance, the "distant reprocessor" model assumes that the 
emitting material in a optically thick shell or ring is initially far enough from the afterglow 
front. Possible conditions of this scenario were discussed in §4.2 (See also Weth et al. 2000; 
Bottcher 2000). On the other hand, the "nearby reprocessor" model requires stellar-like gas 
densities near the GRB site and an ionizing continuum provided by an expanding bubble of 
magnetized plasma breaking through the stellar envelope (Meszaros & Rees 2001), by a pair 
screen scattering a substantial fraction of 7-rays backward (Kumar & Narayan 2002), or by 
continuing emission from the GRB source, e.g., a relativistic wind from a millisecond pulsar 
(Rees & Meszaros 2000). One difficulty when considering the "nearby reprocessor" model 
under the absence of significant line emission arises from our basic ignorance of the continuing 
ionizing source close to the progenitor. This is because the afterglow photons emitted by 
the jet at larger radii render disentangling the contribution close to the progenitor an almost 
infeasible task (Rees & Meszaros 2000). This effect is problematic for the detection of fine 
emission as well since a "photon curtain" of afterglow continuum could swamp any discrete 
line emission. In GRB 020405, the bright continuum could mask line emission occurring well 
within a radius of emission (<C R) even for dense media (n ^ 10^^ cm~^^). In summary, 
a long-lasting, bright afterglow continuum might hinder the detection of faint discrete line 
in the "nearby reprocessor" model unless favorable conditions that might attenuate the 
afterglow continuum (i.e. off-axis explosions, dark bursts?) are in place at the time of the 
observation. 



5. Conclusions and Future Observations 

Our analysis of various discrete X-ray emission models show that the pure power-law 
spectrum of GRB 020405 is adequately explained by a low-density medium with n ^ 5 x 10^^ 
cm~^, or alternatively by a "photon curtain" effect produced by the afterglow continuum in 
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a dense "nearby reprocessor" scenario. In the future, as the rapidity of locahzation improves, 
high-resolution spectroscopy could play a decisive role in resolving the emission mechanism 
(Paerels et al. 2000). 

Figure 4 shows a simulated Gaussian hne profile for Fe Ka with hne width and flux 
identical to those reported for GRB 991216 (Piro et al. 2000), added to the observed power- 
law spectrum of GRB 020405. A Fc Ka line with similar properties would have been clearly 
detected in our LETGS spectrum. But quite possibly the major accomplishment of the Chan- 
dra LETGS in grating spectroscopy could come in the detection of medium-Z a elements, 
Z < 14, where its higher effective area and sensitivity are superior to other instruments on 
board Chandra. Continuous high-resolution spectroscopy might also provide an avenue to 
detect a change in the emission mechanism as a function of time i.e., from recombination 
to fluorescence. In summary, our results confirm the need for further high-resolution obser- 
vations to improve the significance of previous claims of detections for discrete features in 
GRB afterglows. 

We are grateful to all the members of the Chandra team, especially Ed Kellogg, for 
timely planning, implementation and processing of this observation under less than ideal 
circumstances. We also thank Kaya Mori and Maurice Leutenegger for useful conversations. 
Finally we acknowledge Paul Price and the Caltech GRB team for the prompt localization 
of the optical counterpart. 
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Table 1. Line Emission Upper Limits 
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^Rest-frame equivalent width measurement based on a fixed line width model and contin- 
uum level derived from the power-law fit. 

^Indicates lower hmit of recombination coefficient assuming T ^ 4 x 10^ K at the time of 
the observation. 



-14- 




Fig. 1. — Positioning of the Oth-order image on the front-illuminated S2 chip and dispersed 
Ist-order spectra on the S2 and back-illuminated S3 chips in the LETGS/ACIS-S configura- 
tion. 
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Fig. 2.— X-ray light curve of GRB 020405 from the Oth-order LETGS image. The hne 
shows the best-fitted power-law function describing the decay. 
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Fig. 3. — The LETGS/ACIS-S Ist-order coadded spectrum for the whole observation and 
the best-fitted absorbed power law model. 
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Fig. 4. — Simulated Fe Ka line model overlaid on the observed spectrum of GRB 020405. 
Such a feature is clearly ruled out by the data. 



